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Abstract. We study the spatial distribution of satellite galaxies by assuming that they follow the dark matter 
distribution. This assumption is supported by semi-analytical studies based on high-resolution numerical simula- 
tions. We find that for a Milky- Way type halo, if only a dozen satellite galaxies are observed, then they can lie 
on a "great" disk with an rms height of about 40kpc. The normal to the plane is roughly isotropic on the sky. 
These results are consistent with the observed properties of the satellite galaxies in the Milky Way. If, however, 
the satellite galaxies follow the distribution of substructure selected by present mass, then great disks similar to 
the one in the Milky Way are rare and difficult to reproduce, in agreement with the conclusion reached by Kroupa 
et al. (2004). 
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1. Introduction 

The Cold Dark Matter (CDM) structure formation model 
has been successful in explaining many observations, in 
particular the large-scale structure of the universe (e.g. 
Peacock et al. 2001; Tegmark et al. 2004) and the cos- 
mic microwave background (e.g., Spergel et al. 2003). 
However, on small scales, it appears to have problems 
in reproducing the rotation curves of dwarf galaxies and 
the satellite properties in Milky- Way type halos (e.g., see 
Silk 2004 for a review) . But it is still a matter of debate 
whether these conflicts are real and whether they can be 
resolved within the CDM framework. 

The CDM model generically predicts the existence of 
substructures (subhalos). In this model, larger structures 
form by merging of smaller structures, and dense cores 
in these small structures often survive the tidal force 
and manifest as substructures. These substructures are 
most clearly seen in high-resolution numerical simulations 
(e.g., Klypin et al. 1999; Moore et al. 1999; Ghigna et 
al. 2000). If all these substructures host luminous galax- 
ies, then the number of satellite galaxies in a Milky- Way 
type halo can reach several hundred, clearly exceeding 
the number of known satellite galaxies in the Milky Way 
(e.g., Kauffmann et al. 1993). But if only a small fraction 
(~ 10%) of the substructures form stars during their evo- 



Send offprint requests to: X. Kang 
Correspondence to: kangx@shao.ac.cn 



lution, then the number of observed satellites can be rec- 
onciled with observations (Kravtsov et al. 2004). Whether 
the internal kinematics of satellite galaxies are consistent 
with the observations is still unclear (Stoehr et al. 2002; 
Kazantzidis et al. 2004). 

Recently, Kroupa et al. (2004) highlighted the fact that 
the observed spatial distribution of satellite galaxies in the 
Milky Way is distributed in what they term as a "great" 
disk; its plane is almost perpendicular to the plane of the 
Galactic disk. These properties are difficult to understand 
if the satellite galaxies follow the quasi-spherical distri- 
bution expected for the substructures. Recently, Gao et 
al. (2004a; see also Springel et al. 2001; Diemand, Moore 
& Stadel 2004; Nagi & Kravtsov 2005) combined semi- 
analytical techniques and high-resolution numerical simu- 
lations to examine the relation between substructures and 
(satellite) galaxies in clusters of galaxies. They showed 
that satellite galaxies follow roughly the same distribution 
as the underlying dark matter while the substructures fol- 
low a much shallower density distribution in the central 
part. In this paper, we show that if the satellite galaxies 
follow the same spatial distribution as dark matter, then 
many of their puzzling spatial properties highlighted by 
Kroupa et al. (2004) are easier to understand. 

The outline of the paper is as follows. In §2 we first 
describe the simulation data we use, and then compare 
the spatial distribution of satellite galaxies with obser- 
vations, concentrating on the statistics used by Kroupa 
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et al. (2004). We discuss our results in §3. Throughout 
this paper, we adopt the "concordance" ACDM cosmol- 
ogy (e.g., Ostriker & Steinhardt 1995; Spergel et al. 2003 
and references therein), with a matter density parame- 
ter r^m.o — 0.3, a cosmologically constant = 0.7, a 
baryon density parameter fib = 0.024/i~^, and we take 
the power-spectrum normalization erg = 0.9. We write the 
Hubble constant as Hq = 100/ikms~^ Mpc~^ and adopt 
h = 0.7. 



2. Method and Results 

The data we use in this paper is the high resolution halo 
simulation of Jing & Suto (2000). Four halos on galactic- 
mass scales (around 4 x lO^^/i'^Af©) were selected from 
a cosmological simulation with a box size of lOO/i^^ Mpc. 
They were then re-simulated using the nested-grid PPPM 
code which was designed to simulate high-resolution halos. 
The force resolution is typically 0.4% of the virial radius, 
and each particle has mass about 6 x lO^/i^^Af© for the 
four galactic halos. We use the SUBFIND routine of Springel 
et al. (2001) to identify the disjoint self-bound subhalos 
within these halos and all subhalos with more than 10 par- 
ticles are included in our analysis. Furthermore, in order 
to check the distribution of satellite galaxies in galactic ha- 
los, we use the galaxy catalog in a cosmological simulation 
of box 25h^^ Mpc constructed by Kang et al. (2004) using 
their semi-analytical model of galaxy formation. The sub- 
galactic satellites in their simulation were resolved with 
more than 10 particles, and each particle has mass of 
7.7 X IO'^H-^Mq. 

It is known that in galaxy clusters the observed galaxy 
distribution follows the underlying dark matter (Gao et al. 
2004a), but it is less clear whether the observed satellite 
galaxies in galactic mass halos follow the same distribu- 
tion as the dark matter. In Fig.^we show the distribution 
of the observed satellites of the Milky Way and the com- 
parisons with our simulation. The dotted line with trian- 
gles are the model satellite galaxies (with My < —12) in 
galactic halos in our simulation. The thick solid line is the 
distribution of dark matter in the high-resolution galactic 
halos. We found that the modeled satellite galaxies fol- 
low the underling dark matter distribution very well. Also 
it can be seen that the observed satellites roughly follow 
the dark matter particles. The Kolmogorov-Smirnov (KS) 
test shows that their distributions arc consistent at a 60% 
confidence level. In contrast, the substructures with more 
than 10 particles clearly differ from the dark matter or the 
observed satellites distribution. In particular, they are not 
significantly centrally concentrated and have a higher den- 
sity than observation in the outer part, as also found by 
others (De Lucia et al. 2004; Gao et al. 2004a; Diemand et 
al. 2004). Note that our subhalo population is incomplete 
as we can only identify those with more than ten particles 
(> 6 X lO^/i^^Af©). Increasing the mass resolution will 
result in more subhalos, but the relative number density 
profile will not change with the resolution as higher res- 
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Fig. 1. The radial distribution of satellite galaxies. The 
solid circles are the distribution of the 11 inner satellites 
observed in the Milky Way (Mateo 1998; Kroupa et al. 
2004); the Poisson error is indicated for each data point. 
The thick solid line is the dark matter profile while the 
thin solid line is for the substructures. The dotted line with 
triangles is for the satellite galaxies (bright than My = 
— 12) in galactic halos predicted by the semi-analytical 
model of Kang et al. (2004), and the dashed line is for 
the sample selected by mass at accretion (see text). Here 
the radius is normalized to the virial radius and the virial 
radius of the Milky Way is taken to be 250 kpc. 

olution primarily allows the measurement of the density 
profile near the centre (Diemand et al. 2004). 

If all the subhalos form stars then the predicted spa- 
tial distribution of the satellites will be in conflict with 
the observed distribution. Furthermore, many models of 
galaxy formation will over-predict the number of dwarf 
satellite galaxies by an order of magnitude (Kauffmann et 
al. 1993; Bullock et al. 2000; Somerville 2002; Benson et 
al. 2002, Kang et al. 2004). However, if there is a bias 
between the satellite galaxies and the subhalo popula- 
tion, then this over-abundance problem can be alleviated 
(Taylor et al. 2003; Kravtsov et al. 2004). For example, 
using high-resolution N-body simulations, Kravtsov et al. 
(2004) showed that if only a fraction 10%) of the subha- 
los with mass < 10^ — 10^ Af© host luminous galaxies, then 
the observed number of galactic satellites and their spatial 
distribution can be reproduced. They found that these lu- 
minous satellites are descendants of accreted massive halos 
(> lO^M©) at z > 2. These massive systems with tem- 
perature larger than lO^'K can have efhcient gas cooling 
and star formation during their evolutions. Motivated by 
this study, we check whether satellites (in a Milky- Way 
type halo) selected in this way have a similar distribution 
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as the observed satellites. There are about 100 massive 
halos (> 10^h~^ Mq) selected at accretion with redshift 
z > 2. The resolution of our re-simulated halos is not high 
enough for us to fully resolve the dwarf galaxies with mass 
around 10^ Mq. So we, instead, use the dynamical friction 
timescale given by Navarro et al. (1995, see their eq. 2) to 
determine which of the 100 most massive halos can survive 
in the final galactic halo, and we find that 30 of them sur- 
vived. In order to obtain their final positions, we identify 
the most-bound particles of these halos at accretion and 
the final positions of the survied halos at z = are tagged 
by these most-bound particles. This is reasonable as high- 
resolution N-body simulations have shown that the inner 
part of the halo will keep intact during the evolution (e.g. 
Springel et al. 2001). The radial distribution of these sur- 
vived halos is plotted in Figure 1 with the dashed line. 
Clearly the satellite galaxies selected by mass at accretion 
have a similar distribution with the observed satellites in 
the Milky Way. The KS test shows that it also agrees 
with the radial distribution of the dark matter particles 
at a 20% confidence level. We also check if the satellites 
selected in this way have a similar shape with the dark 
matter particles. The number of satellites (30) is not large 
enough to define their shape accurately. Nevertheless, we 
fit their distribution with a triaxial distribution and char- 
acterize their shape as the ratio between the minor and 
major axes. We then use Monte Carlo simulation to pro- 
duce a distribution of such ratios by selecting the satellites 
from the underlying dark matter particles. We found that 
the ratio of the satellites selected by mass at accretion lies 
at the 75% percentile of the distribution. This agreement 
provides the justification to select the satellites galaxies 
from the dark matter particles in the re-simulated galactic 
halos. We also consider the case when the satellite galaxies 
follow the distribution of the subhalos (> 6 x 10^/i^^Mq), 
but we will show that their distribution is not consistent 
with the observations. 

To examine the spatial distribution of satellites in our 
simulated halos, we use Monte Carlo simulations to select 
the satellites assuming they follow either the dark matter 
distribution or the substructure distribution. For a full 
comparison with the results of Kroupa et al. (2004), we 
select 11 random satellites in a sphere with radius of about 
250 kpc from the centre. In order to see how the results 
change with the number of satellites, we also show the 
results by increasing the number of satellites by a factor 
of 2, i.e., 22 satellites within the same sphere. 

Following Kroupa et al. (2004), we fit the selected 
satellites with a plane by minimizing the rms of the height. 
The thickness of the plane is the minimum of the rms of 
the height A. The ratio of A to the maximum distance 
(-Rcut) to the satellites is used to characterize the thickness 
of the plane. In Fig.|21we show the distribution of the char- 
acteristic thickness A/i?cut of the fitted plane from a large 
number of Monte Carlo realizations. The four panels are 
for the four different halos. In these plots the histograms 
without symbols are the distributions of the satellites se- 
lected from the dark matter particles and those with sym- 
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Fig. 2. The distribution of the rms of the scale height 
from a plane defined by satellite galaxies obtained through 
many Monte Carlo realizations. The four panels are for 
four different galactic-sized halos. The histograms with- 
out symbols show that expected if the satellite galaxies 
follow the dark matter distribution while the histograms 
with symbols show those if satellites follow the subhalo 
distribution. The thick solid histograms are for 11 satel- 
lite galaxies, while the thin solid histograms are for 22 
satellite galaxies. The thick horizontal bar in the top left 
panel shows the range of A/i?cut for the observed satel- 
lites when different numbers of satellite galaxies are used 
(see Table 1 in Kroupa et al. 2004). 



bols are for the substructures. We find that for dark mat- 
ter samples, the thickness of the planes are smaller than 
that in the samples of substructures. The average thick- 
ness of the planes peaks at A/i?cut ~ 0.17. It can also be 
seen that increasing the number of satellites also increases 
the thickness of the planes. Kroupa et al. (2004) showed 
that for the inner 11 satellites the characteristic thick- 
ness of the plane is A/i?cut « 0.1, but this A/i?cut value 
changes when one adopts a different number of satellites. 
For example, for the 9 innermost satellites, the character- 
istic thickness of the plane is 0.17. The thick horizontal 
bar in the top left panel of Fig. [21 indicates the range of 
A/i?cut if one adopts 5-16 satellite galaxies in the Milky 
Way. Clearly it is difficult to produce a value as small as 
the observed one if the satellite galaxies follow the sub- 
structure distribution. In contrast, if the satellite galaxies 
follow the dark matter distribution, it is not difficult to 
have the observed low value of A/i?cut- It is easy to show 
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that for a given spherical density profile p oc r the 
thickness of the plane as defined here is given by 

A I 3-n 

t: — = \ / —, r, n < 3. (1) 

i?cut Y 3(5 - n) ' ^ ^ 

The above equation is valid when all the particles in the 
sphere are selected. Notice also that when the power-law 
index n — > 3, the thickness goes to zero. This is because 
the mass and hence the particle number diverge at the ori- 
gin when n — > 3. The plane determined by these (infinite) 
number of particles around the origin will have formally 
zero thickness. 

^From Fig. n we see that the substructures have a shal- 
lower slope (i.e., a smaller value of n) at the centre, so most 
of them are in the outer part of the halo, which makes the 
fitted plane thicker. The observed satellites, in contrast, 
follow a steeper slope than the substructures, so satellites 
are more likely in the inner part of the halo and hence 
make the plane more prominent. In fact the thickness also 
depends on the number of particles selected. Clearly a 
plane with only 3 particles will have zero thickness while 
a plane with infinite number of particles will have a thick- 
ness given by eq. (1). This is why a larger characteristic 
thickness is found when we increase the satellite number 
from 11 to 22 in Fig.|21 

Kroupa et al. (2004) points out that the fitted plane 
of the Milky Way satellites is almost perpendicular to the 
Galactic disk; the normals to these two planes have an an- 
gle of approximately 75°. Unfortunately it is not clear how 
we should define the disk plane in our simulated galaxies. 
There are at least two ways of defining the disk plane. One 
is simply to identify the total angular momentum vector 
as the normal to the disk plane. The second way is that 
we can fit the density profile of the halos with a triax- 
ial distribution, obtain the three principal vectors (a, b 
and c, a < 6 < c), and identify the disk plane as that 
defined by the major and medium axes. These two defini- 
tions would be the same if the minor axis is parallel to the 
angular momentum vector. While there is a statistical cor- 
relation between these two, the scatter is quite large (e.g., 
Faltenbacher et al. 2002). In Fig.|21wc show the angle (/3) 
between the total angular momentum vector and the mi- 
nor axis. Note that we use all the particles enclosed within 
the virial radius to obtain the total angular momentum 
vector and the triaxial density profile. As one can see, the 
alignment for three halos is within 40°, but for the halo in 
the lower left panel, (3 — 78°, i.e., the minor axis is almost 
perpendicular to the total angular momentum vector. The 
identification of the disk plane is further complicated by 
the fact that it is not clear how the angular momentum 
of the (baryonic) disk is related to the total angular mo- 
mentum of the dark matter in N-body/hydro-dynamical 
simulations. Chen et al. (2003) showed that there is a large 
scatter between these two. 

In Fig. O we show the distribution of the angle (9) 
between the normal vector to the fitted plane and the mi- 
nor axis a, assuming that the disk plane is the same as 
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Fig. 3. The distribution of the rms of the angle (in de- 
grees) between the plane defined by satellite galaxies and 
the plane defined by the major and medium axes in the 
triaxial halo model (see text). The dashed line is the pre- 
diction if the angle between the satellite plane and the 
minor axis of the halo is randomly distributed on the sky. 
The other line symbols are the same as in Fig. [3 In each 
panel we label the ratio of the minor, medium, maximum 
axes. In the parentheses we indicate the angle between 
the angular momentum vector and the minor axis of the 
triaxial model. 
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Fig. 4. A scatter plot of 6 vs A/i?cut for the satellites 
selected from the dark matter particles. 6 is the same as 
that in FigEl 



X. Kang et al.: Are Great Disks Defined by Satellite Galaxies in Milky-Way Type Halos Rare in ACDM? 



5 



the plane defined by the major and medium axes. In each 
panel the ratio of a:b:c is also labeled. For two randomly 
oriented normal vectors, the probability distribution for 
the angle between them is oc sin 0, which is shown as the 
dashed line in Fig. |2| Compared with the random distri- 
bution, the satellites are preferentially located in a plane 
defined by the major and medium principal axes. If we 
select only 11 satellite galaxies, then the probabilities for 
the four halos to have an angle larger than 60° (80°) are 
between 28-41% (9%-14%), respectively. So the probabil- 
ity for the fitted satellite plane to have a large angle with 
the galaxy disk is not negligible. On the other hand, if we 
identify the galaxy disk as that given by the angular mo- 
mentum vector, these corresponding probabilities become 
even larger, 28% - 70% (9%-26%) for > 60° {9 > 80°). 

The above results show that there is some probability 
for the selected satellite galaxies to lie in a plane and some 
probability for the plane to be perpendicular to the disk. 
But it is important to examine the joint distribution of 
A/i?cut and the angle 9. Fig. ^ shows the scatter plot be- 
tween A/i?cut and the angle 9 for the satellites from dark 
matter distribution (the scatter plot for substructures is 
similar, but of course with larger values of A/_Rcut)- As we 
can see, there is no correlation between them, for a given 
thickness of the plane, it has roughly the same probabil- 
ity to have large and small values. The probability for 
having > 60° (80°) and A/i?cut < 0.17 ranges from 10% 
to 17% (from 3% to 6%) for the four different haloes. This 
means that it is not rare for the observed plane to have 
large inclinations with the disk. 

Kroupa et al. (2004) have analyzed the cumulative dis- 
tribution of cosw for the Milky Way satellites, where oj is 
the angle between the normal to the great disk and the 
vector connecting the point on the great disk that is clos- 
est to the Galactic centre to the position of the satellite. 
As is designed, the distribution of cosw is closely related 
to the thickness of the great disk. They compared this dis- 
tribution with that of 10^ points derived from a spherical 
power-law density distribution r~". Using the KS test, 
they found the likelihood, Pks, that the observed Milky 
Way distribution is derived from the r~" model is only 
0.005 for < rt < 2.3. Taken at face value, this seriously 
challenges the assumption that the Milky Way satellites 
are consistent with an isotropic parent distribution. 

We have repeated their analysis for the Milky Way, 
and found that Pks is ^ 0.01 if the distribution of | costal 
is used, consistent with Kroupa et al. (2004), who used 
the same quantity although they paper appears to imply 
a different quantity, cosw, was used (Kroupa 2005, pri- 
vate communication). But there is no obvious reason why 
one should use |cosw| instead of cos uj. The distribution 
of cosw should equally well describe the thickness of the 
great disk. When we use coscj, we found that Pks increases 
dramatically to ^ 20% for the Milky Way satellites, which 
implies that the great disk is compatible with the hypoth- 
esis that the satellites are drawn from an isotropic parent 
population. 



We randomly select 11 dark matter particles from the 
first galactic halo, and repeat the above KS test for the 
simulation satellites. We have made many realizations, 
and found that there is a 32% (8.8%) probability that 
the Pks of the simulated satellites is smaller than that 
of the observed cos a; (|cosw|) distribution for the Milky 
Way. About 40% (19%) of the simulated populations have 
Pks < 0.1 for the cosuj ( |cosw|) distributions. These re- 
sults have two important implications. First, the KS prob- 
ability using I coswl (or cosw) should not be directly inter- 
preted as the probability whether the satellites are consis- 
tent with being drawn from a spherical density distri- 
bution, otherwise the probability may be under- estimated. 
The reason is that the great disk is determined using the 
positions of all the satellite galaxies, as a result the cosw 
values based on the great disk are no longer independent 
of each other. This makes the KS probability misleading 
as it assumes that the cos a; values are independent data 
points (e.g., Lupton 1993). Here is an extreme example 
that illustrates this point clearly. If only 3 simulated satel- 
lites are selected from a spherical r~" distribution, then 
cosw is zero for all the three satellites. One finds that 
the KS probability defined above is very small 0.13%). 
This low probability is clearly misleading as this popula- 
tion is drawn from an isotropic power-law density distribu- 
tion. Second, our results clearly show that there is a 32% 
(or 8.8%) probabihty to reproduce the observed cosw (or 
costal) distribution in the ACDM halos, thus it is not so 
rare to expect great disks similar to the one in the Milky 
Way in the ACDM model. 



3. Discussion 

The presence of dark matter substructures is a generic 
prediction of the ACDM model. However, as pointed out 
by many authors (e.g., Springel et al. 2001; Gao et al. 
2004a,2004b; Nagi & Kravtsov 2005), the correspondence 
between dark matter substructures and luminous satellite 
galaxies is not simple. This is because the stellar mass 
of a galaxy may be primarily determined when it is first 
formed. In contrast, the dark matter substructure may 
evolve quite differently due to dynamical friction and other 
processes at later times. In fact, in the study of Gao et 
al. (2004a) , they find ^ 50% of satellite galaxies have no 
corresponding dark matter subhalos. 

If the satellite galaxy distribution follows the dark 
matter distribution, our Monte Carlo simulations demon- 
strate that the spatial distribution of satellite galaxies can 
be better accommodated within the ACDM cosmogony. In 
particular the fact that the Milky Way satellite galaxies 
are distributed in a great disk with its plane almost per- 
pendicular to the stellar disk is not as rare as one naively 
expects. In particular, for the halo shown in the bottom 
left panel of Fig. [SJ as the angular momentum vector of 
the halo is almost perpendicular to minor axis, it has the 
largest probability, > 26% (or 14%) for the fitted plane to 
have an angle larger than 80° with the disk if we identify 
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the disk plane as that given by the total angular momen- 
tum vector (or the minor axis in the triaxial model). 

In contrast, if the satellite galaxies follow the distribu- 
tion of the substructures selected by the present-day mass 
then we find that the satellite galaxy properties will be 
difficult to match because the substructures have a flat- 
ter slope with n ^ 0. But if only ~ 10% percent of the 
substructures host luminous galaxies, the observed distri- 
bution of the satellite galaxies can be recovered, as also 
shown by Kravtsov et al. (2004). It remains, however, 
somewhat puzzling why the observed "great " disk is al- 
most perpendicular to the disk. As we have shown at the 
end of §2, the probability of having a thin great disk per- 
pendicular to the Galactic plane is a few percent to ~ 17%. 
In fact, in our simulations, we find the the fitted plane of 
satellite galaxies have a larger probability to lie in the 
disk defined by the major and medium axes in the triaxial 
model. If the baryonic disk is in the same plane, then the 
great disk in the Milky Way is most easily explained if 
the Galactic plane has a large inclination with the plane 
defined by the the major and medium axes of the dark 
matter distribution. There is currently no observational 
evidence for this. 

Unfortunately, the comparison between our simula- 
tions and observations is not direct. Theoretically, with 
our coUisionless, dark matter only simulations, it is not 
clear how to relate the dark matter properties (such as the 
angular momentum vector) to baryonic disk properties. 
Furthermore, the dark matter halo and satellite galaxy 
profiles will be affected (perhaps adiabatically, e.g.. Mo, 
Mao & White 1998) by the assembly of the baryonic disk. 
Furthermore, the presence of the disk may force subhalos 
into roughly coplanar orbits (Meza et al. 2004) and makes 
it more difficult to have the satellite great disk perpendic- 
ular to the stellar disk. Observationally, it is not clear all 
the satellite galaxies in the Milky Way have been found. 
In particular, some satellite galaxies in low Galactic lati- 
tudes may have been missed (Willman et al. 2004). The 
addition of any such satellites (if at large distances) will 
make the great disk in the Milky Way much thicker. It ap- 
pears that the issue of satellite galaxies can only be fully 
resolved when more observations and high-resolution hy- 
drodynamical simulations become available. 
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